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The influence of baryons on the chiral phase transition in QCD. 
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Abstract 



A qualitative analysis of the chiral phase transition in QCD with two massless quarks and non-zero 
baryon density is performed. It is assumed that at zero baryonic density, p = 0, the temperature phase 
transition is of the second order. Due to a specific power dependence of baryon masses on the chiral 
condensate the phase transition becomes of the first order at the temperature T = rph(p) for p > 0. 
At temperatures rcont(p) > T > Tph(p) there is a mixed phase consisting of the quark phase (stable) 
and the hadron phase (unstable). At the temperature T = rcont(p) the system experiences a continuous 
transition to the pure chirally symmetric phase. 
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(N, 

^^3 I It is well known, that the chiral symmetry is valid in perturbative quantum chromodynamics (QCD) 

^P ' with massless quarks. It is expected also, that the chiral symmetry takes place in fuU-perturbative and 

nonperturbative QCD at high temperatures, (T > 200 MeV), if heavy quarks (c, 6, t) are ignored. The chiral 
symmetry is strongly violated, however, in hadronic matter, i.e. in QCD at T = and low density. What 
i-^ I is the order of phase transition between two phases of QCD with broken and restored chiral symmetry at 

variation of temperature and density is not completely clear now. There are different opinions about this 
Qh! subject (for a detailed review see Ref. ^12 and references therein). 

^5 ' We would like to consider here the influence of baryon density on the chiral phase transition in hadronic 

matter. Kogan, Kovner and Tekin [5] have suggested the idea, that baryons may initiate the restoration 

of chiral symmetry, if their density is high - when roughly half of the volume is occupied by baryons. The 

k>( \ physical argunrent in favour of this idea comes from the hypothesis (supported by calculation in chiral soliton 

j_j ■ model of nucleon |4| ), that inside the baryon the chiral condensate has the sign opposite to that in vacuum. 

Cd I This hypothesis is not proved. Even more, it is doubtful, that the concept of quark condensate inside the 

nucleon can be formulated in a correct way in quantum theory. But the idea on the strong influence of baryon 

density on the chiral phase transition looks very attractive. For this reason no assumption on the driving 

mechanism of chiral phase transition at zero baryon density will be done in this paper. The problem, under 

consideration is: how the phase transition changes in the presence of baryons. The content of the paper 

closely follows ref. 5. 

We discuss the phase transitions in QCD with two massless quarks, u and d. Many lattice calculations 
[HI El IHl EI indicate, that at zero chemical potential the phase transition is of the second order. It will be 
shown below, that the account of baryon density drastically changes the situation and the transition becomes 
of the first order, and, at high density, the matter is always in the chirally symmetric phase. 

Let us first consider a case of the zero baryonic density and suppose that the phase transition from 
chirality violating phase to the chirality conserving one is of the second order. The second order phase 
transition is, generally, characterized by the order parameter 77. The order parameter is a thermal average of 
some operator which may be chosen in various ways. The physical results are independent on the choice of 
the order parameter. In QCD the quark condensate, 77 — \{0\uu\0)\ — |(0|dd|0)| > 0, may be taken as such 
parameter. In the confinement phase the quark condensate is non-zero while in the deconfinement phase it 
is vanishing. 

The quark condensate has the desired properties: as it was demonstrated in the chiral effective theory |10[ 
111) . 77 decreases with the temperature increasing and an extrapolation of the curve 77 (T) to the higher 

temperatures indicates, that r] vanishes at T = Tc « 180 MeV. Here the superscript "0" indicates that 



the critical temperature is taken at zero baryon density. The same conclusion follows from the lattice 
calculations [SI El ^] , where it was also found that the chiral condensate 77 decreases with increase of the 
chemical potential [131 [Tl| . 

Apply the general theory of the second order phase transitions |15| and consider the thermodynamical 
potential <!>(??) at the temperature T near Tc . Since rj is small in this domain, $(77) may be expanded in ry: 

Hv)^^o + Ia7j^ + ^Bt^\ B>{). (1) 

For a moment we neglect possible derivative terms in the potential. 

The terms, proportional to rj and rf vanish for general reasons jl5| . In QCD with massless quarks the 
absence of r] and rf terms can be proved for any perturbative Feynman diagrams. At small t = T — Tc 
the function A{t) is linear in t: A{t) = at, a > 0. If t < the thermodynamical potential $(r?) is minimal 
at 77 ^ 0, while at i > the chiral condensate vanishes 77 = 0. At small t the t-dependence of the coefficient 
B{t) is inessential and may be neglected. The minimum, fj, of the thermodynamical potential can be found 
from the condition, d^/drj — 0: 

r ^/^^tjB, t<0; 
"^'X 0, t>0. ^^> 

It corresponds to the second order phase transition since the potential is quartic in rj and - if the derivative 
terms are included in the expansion - the correlation length becomes infinite at T = Tc . 

Turn now to the case of the finite, but small baryon density p (by p we mean here the sum of baryon 
and anti-baryon densities). For a moment, consider only one type of baryons, i.e. the nucleon. The 
temperature of the phase transition, Tph, is, in general, dependent on the baryon density, Tph = rph(p), 
with T'ph(p = 0) = Tc At T < Tph(p) the term, proportional to Ep, where E = ^yj^ + m? is the baryon 
energy, must be added to the thermodynamical potential |^. As was shown in |16[I17| the nucleon mass m 
(as well as the masses of other baryons) arises due to the spontaneous violation of the chiral symmetry and 
is approximately proportional to the cubic root of the quark condensate: m — cri^'^, with c = (Stt^)^''^ for 
a nucleon. At small temperatures T the baryon contribution to $ is strongly suppressed by the Boltzmann 
factor e~^''^ and is negligible. Below we assume that the proportionality m ~ j]^'^ is valid in a broad 
temperature interval. Arguments in favor of such an assumption are based on the expectation that the 
baryon masses vanish at T = Tpi^^p) and on the dimensional grounds. Near the phase transition point 
E = ^ jp- + rv?' Ki p + c^ rf''^ /{2p). At 7y — > all baryons are accumulating near zero mass and a summation 
over all baryons gives us - instead of eq. QJ - the following: 

^Tl,p) = <^o + \atTf + -^Bi^^ + Cif'^p, (3) 

where C = ^^ cf/(2pi). The term pJ^iPi i^ absorbed into $0 since it is independent on the chiral condensate 
rj. The typical momenta are of the order of the temperature, pi ~ T. Thus, Eq. @ is valid in the region 
T] <^ T^. In the leading approximation the term C can be considered as independent on the temperature at 

Due to the last term in Eq. Q the thermodynamical potential always have a local minimum at 77 = 
since the condensate 77 is always non-negative. At small i < there also exists a local minimum at 7; > 0, 
which is a solution of the equation: 

^^iat + Brj')7^ + ^Cprr'/' = 0. (4) 

At small enough baryon density p, Eq. Q) [visualized in Figure ^^a)] has, in general, two roots, 771 < 770 
and 772 > 770, where rjo — {—at/3B)^'^ is the minimum of the ffi'st term in the right-hand side of Eq. Q. 
The calculation of the second derivative d^^/drf shows that the second root 772 (if exists) corresponds to 
minimum of $(77) and, therefore, is a local minimum of $. The point 7/ = 771 corresponds to a local maximum 
of the thermodynamical potential since at this point the second derivative is always non-positive. 

The thermodynamical potential $(7/,^) at (fixed) non-zero baryon density p has the form plotted in 
Figure ^b). At low enough temperatures (curve Ti) the potential has a global minimum at 77 > and 
system resides in the chirally broken (hadron) phase. As temperature increases the minima at 77 = and 
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Figure 1: Figure 1: (a) Graphical representation of Eq. |0J: "I" is the first term and "11" the second term (with 
the opposite sign) in the r.h.s. of the equation, (b) The thermodynamic potential Q vs. the chiral condensate 
at a fixed baryon density p > 0. At low enough temperatures, T — Ti, the system resides in the chirally 
broken (hadron) phase. The first order phase transition to the quark phase takes place at Tph = T2 > Ti. 
At somewhat higher temperatures, T3 > Tph the system is in a mixed state. The temperature T4 = Tcont 
corresponds to a continuous transition to the pure quark phase, in which the thermodynamic potential has 
the form T5. 



at 77 = ^2 > becomes of equal height (curve T2 = Tph)- At this point the first order phase transition to 
the quark phase takes place. At somewhat higher temperatures, T ^ T^ > Tph, the rj > minimum of 
the potential still exist but $(77 — 0) < $(772)- This is a mixed phase, in which the bubbles of the hadron 
phase may still exist. However, as temperature increases further, the second minimum disappears (curve 
T4 = Tcont)- This temperature corresponds to a continuous transition to the pure quark phase, in which the 
thermodynamic potential has the form T5. 

Let us calculate the temperature of the phase transition, rph(p), at non-zero baryon density p. The 
transition corresponds to the curve T2 in Figuremb), which is defined by the equation $(772, p) = *&(?/ = 0, p), 
where 772 is the second root of Eq. Q as discussed above. The solution is 
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and the second minimum of the thermodynamic potential is at 772 = [4a (Ti'^' ~ Tp\^{p)) / {5 B)] 

At a temperature slightly higher than Tph(p) the potential is minimal at 77 — 0, but it has also an unstable 
minimum at some 77 > 0. The existence of metastable state is also a common feature of the first order phase 
transition {e.g., the overheated liquid in case of liquid-gas system). With a further increase of the density 
p (at a given temperature) the intersection of the two curves in Figure PT^a) disappears and the two curves 
only touch one another at one point rj ~ fJ4. At this temperature a continuous transition (crossover) takes 
place. The corresponding potential has the characteristic form denoted as T4 in Figure^a). The temperature 
T4 = Tcont is defined by the condition that the first iQJ and the second derivatives of Eq. Q vanish: 



Tcont(p)=TW 



2Cp 



3/5 



2/5 



(6) 



and the value of the chiral condensate, where the second local minimum of the potential disappears is given 

(0) 1/2 

by T]4 = [2a(rcont(p) — Tc )/{5B)] . At temperatures T > Tcont{p) the potential has only one minimum 

and the matter is in the state with the restored chiral symmetry. Thus, in QCD with massless quarks the type 

of phase transition with the restoration of the chiral symmetry strongly depends on the value of baryonic 



density p. At a fixed temperature, T < Tc , the phase transition happens at a certain critical density, 
Pph- According to Eq. |(SJ the critical density has a kind of a "universal" dependence on the temperature, 
Pph{T) oc [Tc — r]^", the power of which docs not depend on the parameters of the thermodynamic 
potential, a and B. 
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Figure 2: Figure 2: The qualitative phase diagram at finite baryon density and temperature based on the 
analysis (a) without and (b) with indication of the approximate 2-nd order transition domain. 

The expected phase diagram is shown qualitatively in Figure |2Ia). This diagram does not contain an 
end-point which was found in lattice simulations of the QCD with a finite chemical potential ^lE]- We 
expect that this happens because in our approach a possible influence of the confinement on the order of 
the chiral restoration transition was ignored. Intuitively, it seems that at low baryon densities such influence 
is absent indeed: the deconfinemcnt phenomenon refers to the large quark-anti-quark separations while the 
restoration of the of the chiral symmetry appears due to fluctuations of the gluonic fields in the vicinity of 
the quark. However, the confinement phenomenon dictates the value of the baryon size which can not be 
ignored at high baryon densities, when the baryons are overlapping. If the melting of the baryons happens in 
the hadron phase depicted in Figure |2I a), then at high enough density the nature of the transition could be 
changed. This may give rise in appearance of the end-point observed in Ref. 18.. 19; . The domain where the 
inequality \at\ 3> Cprf^^^ p ^ is fulfilled, has specific features. In this domain the phase transition looks 
like a smeared second order phase transition: the specific heat has (approximately) a discontinuity at the 



■ri°)y 



-1/2 



at T - T, 



(0) 



0. 



phase transition point, ACp = a?Tc/B. The correlation length increases as (T ■ 
The latter arises if we include the derivative terms in the effective thermodynamical potential. The phase 
diagram with this domain indicated may look as it is shown in Figure [2Ib). Note that the applicability of 
our considerations is limited to the region |T — Ti \/Tc <C 1 and low baryon densities. 

In the real QCD the massive heavy quarks (the quarks c, 6, t) do not infiuence on this conclusion, since 
their concentration in the vicinity oiT k. Tc ^ 200 MeV is small. However, the strange quarks, the mass 
of which rus ~ 150 MeV is just of order of expected Tc , may change the situation. This problem deserves 
further investigation. 
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